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ABSTRACT 

Context. IGR J19294+1816 was discovered by INTEGRAL in 2009 during a bright X-ray outburst and was classified as a possible 
Be X-ray binary or supergiant fast X-ray transient. 

Aims. On 2010 October 28, the source displayed a second X-ray outburst and a two months-long monitoring with Swift was carried 
out to follow the evolution of the source X-ray flux during the event. 

Methods. We report on the INTEGRAL and Swift observations of the second X-ray outburst observed from IGR J19294+1816. 
Results. We detected pulsations in the X-ray emission from the source at ~I2.5 s up to 50 keV. The source X-ray flux decreased 
smoothly during the two months of observation, displaying only marginal spectral changes. Owing to the relatively rapid decay of 
the source X-ray flux, no significant variations of the source spin period across the event could be measured. This prevented a firm 
confirmation of the previously suggested orbital period of the source at 117 d. This periodicity was also searched for in archival 
Swift [BAT data. We detected a marginally significant peak in the periodogram and determined the best period at 116.2+0.6 days 
(estimated chance probability of a spurious detection 1%). 

Conclusions. The smooth decline of the source X-ray flux across the two months of observations after the onset of the second outburst, 
together with its relatively low value of the spin period and the absence of remarkable changes in the spectral parameters (i.e., the 
absorption column density), suggests that IGR J19294+I8I6 is most likely another member of the Be X-ray binaries discovered by 
INTEGRAL and not a supergiant fast X-ray transient. 

Key words. X-rays: binaries - stars: individual IGR J19294+I816 - stars: neutron - X-rays: stars 



1. Introduction 

IGR J 19294-1- 18 16 was discovered undergoing a bright X-ray 
outburst for the first time by INTEGRAL on 2009 March 27 
dTurleretaLl 12009). On this occasion, IGR J19294H-1816 was 
detected by the IBIS/ISGRI telescope at a significance level of 
8.3 cr in the 20-40 keV energy band. The source was not detected 
in the 40-80 keV energy band. The IBlS/ISGRl spectrum (18-50 
keV) could be reasonably well fitted with a power-law model. 
The estimated photon index was r=4.0±0.7 and the correspond- 
ing flux 14±2 mCrab (i.e., 1.1x10"'" erg cm"^ s"'). 
,—1 The analysis of archival Swift observations permitted obtain- 
i ing a refined position for the source at a'j2ooo=19''29™55'!9 and 
■ 5j20oo=18°18'38'.'4 (associated uncertainty of 3.5'.' at 90% c.l.), 
and thus identify the IR counterpart to IGRJ19294-I-1816 as th e 
object 2MASSJ19295591 + 1818382 (Rodrig uez etaPbOOgbh . 
_ These data also showe d evidence of pulsation s at ~12.5 s, later 
confirmed with RXTE dStrohmaver et al.ll2009l) . 

The historical Swift [BAlT lightcurve of the source revealed 
the presence of a modulation with a period of 117 d, which 
was tentatively interpreted as the orbital p eriod of the system 
dKrimm et alJl2009HCorbet & Krimmll2009n . 

A summary of all these results and an in-depth inspec- 
tion of all available INTEGRAL observatio ns performed in 
the dire ction of the source was provided by [Rodriguez et al.l 
(l2009ah . The authors reported on the discovery of relatively 
short flares (few thousands of seconds) from this source with 
INTEGRAL, and concluded that, even though the orbital and 
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spin period of IGR J19294-I-1816 would nicely place the source 
in the region of the Pspin-^'orb diagram populated by the Be X- 
ray binaries, a possible supe rgiant fast X-ray (SFXTs, see e.g., 
IWalter & Zurita Herasll2007l) nature of this source could not be 
ruled out. 

IGR J19294+1816 was detected undergoin g a renewed ac- 
tivit y with INTE GRAL on 2010 October 28 ( Bozzo et al.ll2010l: 
Jen ke et al.ll2010h . Motivated by the previous findings, we trig- 
gered a monitoring program of this source with Swift /XRT. This 
permitted us for the first time to follow in detail the evolu- 
tion of the source X-ray flux for about two months after the 
new outburst. In this paper, we report on the results of the 
INTEGRAL/lSGRl and Swift /XRT obsei-vations earned out 
during the outburst in 2010 and reanalyze all previously avail- 
able observations performed with these instruments to investi- 
gate the real nature of the source. In Sect. |2]we describe the data 
analysis technique and the results. Our discussion and conclu- 
sions are summarized in Sect. |3] 



2. Data analysis 

2.1. INTEGRAL 

INTEGRAL observations are commonly divided into "science 
windows" (SCWs), i.e., pointings with typical durations of ~2- 
3 ks. We considered all avai l able SCWs for flie IBIS/ISGRI (17- 
80 keV, iLebrun et al.ll2003l: lUbertini et an i2003) telescope that 
were performed in the direction of IGR J19294-I-1816 close to 
the period in which a renewed activity from this source was 
discovered (see Sect. [TJ. These comprised observations of the 
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Fig.l. Mosaic of the IBIS/ISGRI FOV around 
IGRJ19294+1816 realized by using all available SCWs in 
the direction of the source during satellite revolution 982 
(20-50 keV energy range). 
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Fig. 2. IBIS/ISGRI spectrum of IGR J 19294+ 18 16 extracted 
from the observations in Rev. 982. The best fit is obtained with a 
simple power-law model and is shown in the figure with a con- 
tinuous line. The residuals from this fit are reported in the bottom 
panel. 



region around GRS1915 + 105 in satellite revolutions 980, 982, 
985, 986, and 988. A complete log of the observations is pro- 
vided in Table [1] The source was always outside the field of 
view (FOV) of the lower ene rgy telescope JEM-X on-board 
INTEGRAL dLund et al.ll2"003l) . All INTEGRAL data were an- 
alyzed u^ingj[ersi£n_9^ofJhe OSA software distributed by the 
ISDC (ICourvoisier et al.ll2003h . 

IGR J19294+1816 was detected in rev. 982 at a significance 
level of 27.3 cr in the 20-50 keV energy_band (this energy band 
was chosen to maximize the S/N, see fBozzo et a l. 2010). An im- 
age of the IBIS/ISGRI FOV around IGR J19294+1816 is shown 
in Fig.[T] Given the relatively high detection significance, we re- 
binned the IBIS/ISGRI response matrixO to extract a spectrum 
with 22 energy bins and look for possible spectral features. A 
reasonably good fit to this spectrum was obtained using a sim- 
ple power-law model (x^^Jd.o.f. = 1.4/18). We measured a power- 
law photon index of r=4.9+0.3 and estimated a source flux of 
(5.4±0.3)xl0-i° erg/cm^/s (20-50 keV). The IBIS/ISGRI spec- 
trum of IGR J 19294+ 18 16 is shown in Fig. |2] together with the 
best-fit model. No convincing spectral features emerged from the 
residuals of this fit. 

We also analyzed the closest INTEGRAL observations per- 
formed in the direction of the source before the discovery of 
its renewed activity. IGR J19294+1816 was in the IBIS/ISGRI 
FOV from 2010 October 22 at 23:04 to 2010 October 23 at 
02:21 (UTC; satellite revolution 980). In this case the source 
was not detected. We estimated a 3 cr upper limit on its X-ray 
flux of 5.1 mCrab (20-50 keV energy band; effective exposure 
time 22 ks). This correspond^ to 5x10"" erg cm"^ s~^ 

A similar analysis was also performed on the IBIS/ISGRI 
mosaic realized by using the available INTEGRAL observa- 
tions including IGR J19294+1816 in their FOV carried out im- 
mediately after the discovery of the renewed activity (satellite 
Rev. 985, 986, 988). During this period (from 2010 November 



' See http://isdcu l3.unige.ch/Soft/download/osa/osa-do^ osa_doc- 
9.0/osa_um_ibis-9.2/. 

^ The conversion between the source count-rate and mCrab is car- 
ried out by using the IBIS/ISGRI observation of the Crab in Rev. 967. 
Here the count-rate of the Crab in the 20-50 keV energy band was 
168.86+0.09 and the corresponding flux 1.0x10"* erg cm"^ s"' 
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Fig. 3. Power spectrum of the IBIS/ISGRI data of 
IGRJ19294+1816 in Rev. 982. Pulsations at a frequency 
of 0.080312(6) Hz (eiTor at 90% c.l. on the last digit) are cleai'ly 
detected. 



6 at 17:46 to 2010 November 18 at 00:26 UTC) the source was 
not detected and we estimated a 3 cr upper limit on its X-ray 
flux in the 20-50 keV energy band of 3.1 mCrab (i.e. roughly 
3x10"" erg cm"^ s"'; effective exposure time 56 ks). 

Pulsations were searched in the INTEGRAL data of 
rev. 982 by extracting source and background lightcurves with 
time resolutions of 1 s applying the method described in 
Segreto & Ferrigno] (|2007[ all event times were also corrected 
to the solar system barycenter). Pulsations were clearly detected 
in these data (see Fig.|3]l. We estimated the best pulse frequency 
at 0.080312(6) Hz (error at 90% c.l on the last di git) using stan- 
dard phase-shift techniques (Ferrigno et al.l20O7 ^. and measured 
an rms fractional amplitude of 0.84+0.06 (defined according to 
Eq. 14 in lvanderKlislll988l) . The IBIS/ISGRI pulsed profile of 
IGR J 1 9294+ 1 8 1 6 is shown in Fig. S 

In order to better compare these results with the previous 
findings, we also reanalyzed the INTEGRAL data of the pre- 
vious outburst of IGRJ19294+1816 (Rev. 788, see Sect. [B- 
In this case, we found no clear indication of pulsations in the 
IBIS/ISGRI periodogram (the source X-ray flux was a factor 



E. Bozzo et al.: IGR J19294+1816 in outburst 



3 



Rev. 982 



^ o - 

O 
O 

5 

o ^ - 

s 



D.5 1 1.5 

phase 

Rev. 788 



> 

o ~_ 

O 

S - 
O 

s 



0.5 1 1.5 2 

phase 

Fig. 4. IBIS/ISGRI lightcurve of IGR J 19294+ 18 16 in rev. 982 
(top panel) and rev. 788 (bottom panel) folded at the best de- 
termined spin frequency. In both cases we used data in the 20- 
50 keV energy range. 



Table 1. INTEGRAL observations log of IGR J19294+1816. 



Rev. 


START TIME 


STOP TIME 


EXP 


Fx 




(MJD) 


(MJD) 


(ks) 


(erg cm"^ s"') 


788 


54917.49852 


54917.73836 


12 


1.1x10-'" 


980" 


55491.96168 


55492.22987 


22 


<5.0xlO-" 


982 


55497.76780 


55498.00752 


3.4 


5.4x10"'° 


985,6,8" 


55506.74066 


55518.01843 


56 


<3.0xlO-" 



NOTE: Fx is the flux in the 20-50 keV energy band. EXP indicates the 
elfective exposure time of each obsereation. " : 3o" upper limit. 



of ~5 lower than that in rev. 982). However, an epoch-folding 
search around the previously know period led to a marginal de- 
tection of pulsations. We determined a best pulse frequency of 
0.08036(2) Hz, and estimated a fractional rms of 0.55+0.12. The 
IBIS/ISGRI lightcurve of this observation folded at the best de- 
termined spin frequency is shown in Fig. |4] 

A log of all the INTEGRAL observations used is reported in 
Table H 

2.2. Swift 

Following the detection of a new outburst from 
IGR J 19294+ 18 16, we requested a monitoring observation 



of the source to follow closely the evolution of its X-ray flux 
down to the quiescent level. A log of all Swift /XRT observations 
is provided in Table |2] 

Swift PiKT data w ere analyzed by using standard procedures 
(iBurrows et al.l l2005l) and the latest calibration files available. 
The XRT data were processed with the xrtpipeline (v.0.12.6); 
filtering and screening criteria were applied by using ftools 
(Heasoft v.6.10). We extracted source and background light 
curves and spectra by selecting event grades of 0-2 and 0-12 for 
the window time and photon counting (PC) mode, respectively. 
Exposure maps were created through the xrtexpomap task, and 
we used the latest spectral redistribution matrices in the heasarc 
caUbration database (v.Ol 1). Ancillary response files, accounting 
for different extraction regions, vignetting and PSF corrections, 
were generated by using the xrtmkarf task. When required, we 
corrected PC data for pile-up, and used the xrtlccorr to account 
for this correction in the background-subtracted light curves. 
Source and background event lists were extracted by using the 
highest possible time resolution for the PC (2.5 s) and the WT 
(0.0018 s) observational modes, and then barycentered by using 
the barycorr tool available within the Heasoft software pack- 
age. 

From the observations in Table |2]that were characterized by 
adequate statistics, we extracted source and background spec- 
tra and rebinned these spectra to have at least 20 photons per 
bin to permit fitting. We derived a mean X-ray flux by fitting 
these spectra with an absorbed power law model (we used Xspec 
V. 12.6.0). The spectra extracted during observations character- 
ized by a lower statistic were instead rebinned to have at least 
five photons per bin and were then fitted by using the C-statistic. 
For the observations in which less than 100 counts were recorded 
from the source, we estimated the source count-rate with sosta 
(xEviAGE V.4.5.1) and then used this count-rate within webpimm^ 
to derive the X-ray flux (we assumed the same spectral model of 
the closest observation for which a proper spectral analysis could 
be carried out; see e.g., Bozzo et al. 2009) . All these results are 
reported in Table |2]and used in Fig.|5]to plot the variation of the 
source X-ray flux during the decay from the outburst. 

Even though the observations were not strictly simultane- 
ous, we performed a broad-band (0.3-50 keV) fit to the source 
spectrum by combining the IBIS/ISGRI data from rev. 982 and 
the Swift piRT data from the observation ID. 00031854001. We 
used first a simple absorbed blackbody or power-law model, and 
introduced a normalization constant to take into account both 
file intercalibration between the Swift /XRT and IBIS/ISGRI in- 
struments and the variability of the source. Both fits gave un- 
acceptable results with x^g^>'^-^- We therefore substituted the 
power-law component with a cut-off power law (wabs*cutoffpl 
in xspec). This model gave an acceptable fit Oif^^d.o.f.^l.l/Vl) 
and we measured a power-law photon index of r=0.4+0.3, an 
absorption column density of A^h=(3.1±0.7)x10^^ cm"^ and a 
cut-off energy of 8.0;^[ ^ keV. The normalization constant turned 
out to be 4.2;^^] 2, which indicates (as expected) that the flux of 
the source during the Swift observation already decreased with 
respect to that measured by INTEGRAL during satellite rev. 982 
(see also Tab.[T]and|2|. The unfolded XRT+ISGRI spectrum, the 
best-fit model and the residuals from this fit are shown in Fig. |6] 

We used the Z^-st atistic technique (iBuccheri et al.l 119831: 
iMarkwardt et al.l l2002l) to search for pulsations in all the 
Swift PCRT observations reported in Table |2] In each of these 
observations the value of the source spin frequency was es- 

^ http://heasarc.nasa.gov/Tools/w3pimms.html 
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Table 2. Swift observations log of IGRJ19294+1816. 



OBSID 


INSTR 


START TIME 
(MJD) 


STOP TIME 
(MJD) 


EXP 

(sec) 


Nii 

(10^2 cm-2) 


r 


fobs 

(erg cm" 


2 s-') 


(C-stat/d.o.f.) 


2007 


00037191001 


XRT/PC 


54443.05208 


54444.94375 


7.9 


4 4+1-1 
^- -0.8 


OKI 


3.2x10- 


11 


0.9/70 


00037191002 


XRT/PC 


54447.01389 


54447.35347 


3.4 


3 7+'-3 
^■'-1.1 


0.9±0.4 


3.2x10- 


11 


0.9/33 


2009 


00031392001 


XRT/PC 


54925.80417 


54926.00000 


2.6 


^■"-3.6 


1 l + '-5 


6.5x10- 


12 


(17.4/19) 


00031392002° 


XRT/PC 


54973.86319 


54973.93889 


2.0 


3.6 (fixed) 


1.1 (fixed) 


2.0x10- 


12 




00031392003° 


XRT/PC 


54974.53819 


54974.67153 


1.7 


3.6 (fixed) 


1.1 (fixed) 


2.2x10- 


12 




00031392004* 


XRT/PC 


54978.02153 


54978.16319 


1.2 


3.6 (fixed) 


1.1 (fixed) 


1.3x10- 


12 




2010 


00031854001 


XRT/PC 


55498.62225 


55498.63470 


0.99 


o 7+0.8 
•'■'-0.7 


1.1±0.3 


2.6x10- 


10 


1.8/52 


00031854002 


XRT/WT 
XRT/PC 


55498.69269 


55498.88016 


0.12 
2.6 


4.4+2.3 

3 0+''-^ 


1 5+0-^ 

^■■'-0.6 
1.0±0.2 


2.6x10- 
2.3x10- 


10 
10 


0.8/11 
0.7/72 


00031854003 


XRT/PC 


55500.21293 


55500.22519 


0.95 


10 7+^ -° 


0.7±0.7 


6.2x10- 


11 


(46.8/50) 


00031854004 


XRT/PC 


55501.63127 


55501.64450 


0.97 


7 4+3-' 


^- -0.6 


4.5x10- 


11 


(48.1/48) 


00031854005 


XRT/PC 


55503.63889 


55503.97292 


1.3 


7 0+'* ' 
'■"-3.4 


1 0+"-' 


3.0x10- 


11 


(34.8/41) 


00031854006° 


XRT/PC 


55505.50556 


55505.71944 


1.4 


7.0 (fixed) 


1.0 (fixed) 


8.2x10- 


12 




00031854007 


XRT/PC 


55507.79097 


55507.92986 


1.0 


4 8+2-8 

^■°2.2 


1 5+"-'* 

'■■'-0.4 


2.5x10- 


11 


(43.0/36) 


00031854008 


XRT/PC 


55510.45347 


55510.53889 


2.4 


4 7+4.0 
■ -3.6 


1 1 + 1-" 
'■ -0.9 


1.5x10- 


11 


(17.0/27) 


00031854009° 


XRT/PC 


55513.46458 


55513.54931 


2.2 


4.7 (fixed) 


1.1 (fixed) 


4.3x10- 


12 




00031854010° 


XRT/PC 


55516.56111 


55516.63472 


2.3 


4.7 (fixed) 


1.1 (fixed) 


1.0x10- 


11 




00031854011 


XRT/PC 


55523.03056 


55523.10972 


2.3 


"■^-4.2 




8.1x10- 


12 


(20.0/24) 


00031854012 


XRT/PC 


55530.05625 


55530.26944 


2.8 


4 8+2-3 


1.2±0.7 


9.3x10- 


12 


(28.0/31) 


00031854013 


XRT/PC 


55537.01806 


55537.10139 


2.2 




2.8+2-' 


3.4x10- 


12 


(13.2/11) 


00031854014 


XRT/PC 


55544.10764 


55544.99028 


3.3 


4 4+3.0 
■ -3.2 


1.3!j:? 


6.9x10- 


12 


(19.4/24) 


00031854015 


XRT/PC 


55550.47847 


55550.62917 


3.4 


4 6+2-9 

■ -2.5 


1 9+0.9 
'■ -0.8 


5.6x10- 


12 


(11.8/21) 



NOTE: Spectra extracted from these observations are fitted with an absorbed power law (absorption column density Nh and photon index F). 
fobs is the absorbed flux in the 1-10 keV energy band. EXP indicates the total exposure time of each obsereation 
°: We used sosta and webpimms to estimate this flux. *: 90% c.l. upper Hmit. 



timated accurately by fitting the peak in the Z^-statistic peri- 
odogram with a model comprising a constant and a sync func- 
tion. However, because the points in the Z^-statistic periodogram 
are not independent, the error on the position of the centroid of 
the peak derived in this way could not be considered as a re- 
liable estimate of the uncertainty on the source spin frequency. 
We thus simulated for each observation 100 event files with the 
same exposure and number of events as the real file (we assumed 
a constant flux during each observation). A sinusoidal signal at 
the frequency determined from the data was introduced artifi- 
cially and then measured a posteriori by using the Z^-statistic 
(we investigated in each case a total of 1000 frequencies in the 
range 0.078-0.082 Hz). The frequencies with a detection signif- 
icance >3cr determined with this method were then averaged to 
derive the value of the best spin frequency for each observation 
and the associated uncertainty at Icr c.l. (that we assumed equal 
to 1 standard deviation from the averaged value). We checked for 
each observation that the frequencies determined from the simu- 
lations with this method were all consistent within 1 standard de- 
viation with the best value found by using the Z^-statistic on the 
real event file. The results of this procedure are shown in Fig . |7l 
In the bottom panel of this figure, we also report for each obser- 
vation the percentage of the 100 simulated lightcurves that did 
not result in a significant detection of a periodicity (<3cr). The 



observations characterized by a lower source flux or a shorter 
exposure gave, as expected, a higher percentage of failed detec- 
tions. This suggests that the percentage reported in the figure can 
be considered a reliable estimate of the detection significance of 
the periodicity in each observation. 

The results reported in Fig. |7] show that owing to the rel- 
ative large uncertainties that affect most of the measurements, 
no significant variations of the source spin frequency could be 
measured during the Swift follow-up observations. This pre- 
vented any even tentative determination of an orbital solution 
for IGR J19294-I-1816 through the pulse arrival time technique. 

To search for a confirmation of the suggested ~117 d orbital 
period of IGR J 19294-1- 18 16 (see Sect.Q, we therefore adopted 
an alternative method. We obtained the long-term lightcurve 
of the source from the Swift [BAT survey web-pagqj binned 
snapshot-by-snapshot. All photon arrival times were corrected to 
the solar system barycenter with the task earth2sun. Because the 
Z^-statistic technique is particularly suited for un-binned event 
files rather than for lightcurves, we searched for periodicities 
in the BAT data by using the Lomb-Scargle ( LS) periodogram 
(we adopted the fast algorithm described in ^ Press & Rvbickil 
1989, and scanned periods above 31 days with an oversampling 



http://heasarc.gsfc.nasa.gov/docs/swift/results/bs58mon/. 
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Fig. 5. Lightcurves of the outburst decay of IGR J 19294+ 18 16 
during the two events that occurred in 2009 (lower panel) and 
2010 (upper panel). We did not plot the observations carried out 
in 2007 because in that case the beginning of the outburst was 
not caught by any X-ray satellite (see Tables[T]and|2]i. We scaled 
the values on the X-axis to the first INTEGRAL detection of the 
source (first point on the left in both panels). This corresponds to 
55497.8877 MJD for the upper panel and 54917.6184 MJD for 
the lower panel. 
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Fig. 7. Results of the search for pulsations in the Swift /XRT data. 
The upper (middle) panel show the determined frequency (pe- 
riod) of pulsations from the source in each observation as a func- 
tion of time. In the bottom panel we show the percentage of the 
total number of simulated lightcurves (100) for each observation 
that did not result in a statistically significant (>3cr) detection 
of a spin period at the expected frequency. The reference time 
(f = 0) is 55498.62225 MJD. 
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Fig. 6. Unfolded joint Swift /XKY and INTEGRAL /ISGRl spec- 
trum (we used the XRT spectrum extracted from observation 
ID. 00031854001 and the ISGRI spectrum extracted from the 
INTEGRAL rev. 982). The best-fit model is obtained with a 
CUTOFFPL model (see text for details). The residuals from this 
fit are shown in the bottom panel. 



factor of 8). Because the uncertainties affecting the Swift fBA^ 
data change significantly from pointing to pointing because 
of the different exposure time and off-axis position of the 
source with respect to the instrument aim-point, we also imple- 
mented the_weighting scheme of the BAT lightcurve proposed 
by Corb et et al.l (l2007l in our case the correction factor Vs intro- 
duced to account for the source variability turned out to be neg- 
ative and thus was set to zero). The Lomb-Scargle periodogram 
was realized for both the original Swift /BAlT lightcurve, and for 
the 500 lightcurves simulated from it. In performing the simula- 
tion, we assumed for each lightcurve a duration equal to that of 
the true Swift /BAT lightcurve and drew the source count-rate in 



each time bin from a normal distribution centered on the real 
measured count-rate0. We set the corresponding standard de- 
viation equal to the measured error on the count-rate. The re- 
sults of this analysis are reported in Fig. |8] The histogram of 
the periods determined from the simulated data showed a sin- 
gle prominent peak at ~120 days, which thus turned out to be 
the most probable period (the simulation does not assume a pri- 
ori a periodicity). However, we noticed that only 266 out of 
the 500 simulated lightcurves displayed a peak at a period of 
~120 days, and from these we estimated the best orbital period at 
f orb=120+2 days (68% c.l.). Most of the remaining lightcurves 
(209) show a prominent peak at frequencies close to higher har- 
monics of the fundamental frequency. This distribution of the 
power among se veral harmonics would b e expected according 
to the findings o f ICorbet & Krimml (l2009h . who showed that the 
folded lightcurve of the source was strongly not sinusoidal. 

To refine our determination of the orbital period, we also 
looked for periodicities in the Swift /BAT lig htcurve by exploit - 
ing the epoch-folding technique described bv lD'Ai et al.l (1201 lb . 
We folded the light-curve in 16 phase bins by using periods in the 
range 15-600 d. The sampling in period was chosen to be 3 times 
finer than the resolution of the method, i.e. Av-P^/(NAsat) 
(here, P is the period, is the number o f phase bins, and Abat 
is the data span length; iBuccheri et al.l [1983). We then com- 
puted for each of the folded lightcurves the corresponding and 
found the most prominent excess at a period of 1 16.2 + 0.6days. 
The estimated chance probability that this detection would corre- 
spond to a spurious effect is about 1% (see Fig.|9]l. In the bottom 
panel of Fig. |9]we also show the Swift /BAlT lightcurve folded at 
the best determined spin period. The reference epoch was chosen 
at 54877.3957 MJD, and we notice the presence of a prominent 
peak in the source count rate at phase 0.34+0.03. Considering an 



' This assumption is justified for a coded mask instrument, because 
the source signal typically constitutes only a few percent of the large 
number of counts that characterize the background. 
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Fig. 8. Results of the simulations carried out to investigate 
the significance of the period at ~120 days found from the 
Swift /BAlT data. The upper panel shows the value of the period 
for each simulation (the number of the simulation is reported on 
the X-axis). The middle panel reports the probability that the de- 
tection of a period is false in each simulation (i.e., a value of 
0.1 corresponds to a detection significance of 90% c.l.), and the 
lower panel the power of each peak in the corresponding LS di- 
agram. 

orbital period of 1 16.2 + 0. 6 days, the epoch of the peak would 
agree with the prediction of ICorbet & Krimml (l2009h . 

We conclude that further observations are probably needed 
to firmly establish and eventually refine the orbital period of the 
source. 

3. Discussion and conclusions 

We reported on the INTEGRAL and Swift observations car- 
ried out for about two months after the newly discovered out- 
burst of IGRJ19294+1816 on 2010 October 28 (Bozzo et al. 
12010 ). During these observations, IGR J19294H-1816 displayed 
a smoothly decreasing X-ray flux with a decay time-scale simi- 
lar to that observed during the outburst in 2009. Pulsations were 
clearly detected in the INTEGRAL data up to ~50 keV and in 
the Swift PCRT data. Owing to the relatively low X-ray flux of 
the source during most of the XRT observations, the accuracy 
with which we could measure its spin period was far too low 
to determine an orbital solution for IGR J 19294+ 18 16 with the 
pulse arrival times technique. Instead, we found a marginal de- 
tection of this periodicity with both the LS periodogram and 
the epoch-folding technique. We found the best orbital period 
at 1 16.2+0.6 days, with a chance probability of 1% to be a spu- 
rious dete ction. 

As Rodr iguez et al.l (l2009ah remarked, the spin and orbital 
period of IGRJ19294+1816 (if confirmed) would nicely place 
the source in the part of the /'spin-^'oib diagram ([Corbet 1986) 
populated by the so-called Be X-ray binaries. The INTEGRAL 
and Swift observations of the source carried out during the out- 
burst in 2010 provide further support to this hypothesis. The 
X-ray pulsations at ~12.5 s and the high rms fractional am- 
plitude measured up to 50 keV are typical of Be X-ray bina- 
ries and are not usually observed in the SFXTs (the ~6 days 
long outburst shown in Fig. [5] would also be typical for a Be 
X-ray binary; see e.g., IZioIkow ski 2002; Reig 2011, for re- 
cent reviews). Among all the candidate and confirmed sources 
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Fig. 9. Results of the search for periodicity in the SwiftfBA T 
data obtained with the technique described bv iD' Ai et al.l(l201 ih . 
The upper panel shows the computed from the Swift /BAlT 
lightcurve folded at different periods in the range 15-600 days. 
The most prominent peak is at 1 16.2 ± 0.6 days. To estimate the 
si gnificance of this peak, we followed the procedure described 
in lD'Ai et a Ddmi]), and report the histogram of the detrended 
values in the middle panel. We show in the bottom panel the 
Swift /BAlT lightcurve (15-50 keV) of the source folded at the 
best determined orbital period (1 16.2 ± 0.6 days). The reference 
time used is 54877.3957 MJD. Each phase bin in the plot corre- 
sponds to an integration time of three days. 



in the SFXT class (~15 objects), so far onl}Q the intermediate 
SFXT IGR J 18483-03 11 showed pulsations with a period simi- 
lar to fliat of IGRJ 19294+1816 (-21 s detected up to -40 keV, 
ISguera et al.ll2007h . However, the smooth and gradual decline 
of the X-ray flux from IGR J 19294+ 18 16 during two months 
spanned by the Swift observations seems to be different from the 
behavior displayed by the SFXT sources. Indeed, these sources 
usually show a rapid (thousands of seconds) decline of the X- 
ray flux of 3-5 orders of magnitude after the occurrence of an 
outburst. In those cases in which a detailed monitoring could 
be carried out with Swift /XRT after the occurrence of a bright 
event, it was also noticed that a pronounced variability in the 
X-ray flux (a factor ~100), accompanied by significant changes 



* Note that Bozzo et al. (2011, submitted) showed that the NS hosted 
in the SFXT IGR J 184 10-0535 might not be pulsating at -4.7 s as it was 
previously reported. 
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in the spectral parameters (e.g., the absorption column density), 
characterize the source emission during the entire orbit (see e.g., 
[Bozzo et al. 2008, 2009; Sidoli et al. 2009). The slow decline of 
the X-ray flux from IGR J 19294+ 18 16, combined with the fairly 
small changes in the spectral parameters across the entire Swift 
monitoring, suggest that the behavior of IGR J19294+1816 is 
closer to that of the Be X-ray binaries. We thus conclude that 
IGR J 19294-1- 18 16 is most likely a Be X-ray binary and not an 
SFXT. 
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